ctf 



-a 

o 
o 



> 

o 



Condensed Matter Physics, 201 1 , Vol. 14, No 2, 23701 : 1-01 QSffiQHnMlJD 

DOI: 10.5488/CMR1 4.23701 [MMWHi! 

http://www.icmp.lviv.ua/journal Lr ^ LnJ J ^- )LI - J ^ J 



A comparative study for structural and electronic 
properties of single-crystal ScN 

R. MohammacPEl, §. KatirciogliPl 

1 Palestine Technical University, Applied Science College, WestBank, Palestine 

2 Middle East Technical University, Physics Department, 06530 Ankara, Turkey 

o 

{vq ■ Received November 4, 201 0, in final form February 3, 201 1 

A comparative study by FP-LAPW calculations based on DFT within LDA, PBE-GGA, EV ox -PW co -GGA, and 
EV ox -GGA-LDA co schemes is introduced for the structural and electronic properties of ScN in RS, ZB, WZ, 
and CsCI phases. According to all approximations used in this work, the RS phase is the stable ground state 
structure and makes a transition to CsCI phase at high transition pressure. While PBE-GGA and EV ox -PW co - 
GGA's have provided better structural features such as equilibrium lattice constant and bulk modulus, only 
EVex-PW co -GGA and EV cx -GGA-LDA co 's have given the non zero, positive indirect energy gap for RS-ScN, 
~~ ' comparable with the experimental ones. The indirect band gap of ScN in RS phase is enlarged to the corre- 

sponding measured value by EV cx -PW co -GGA+U SIC calculations in which the Coulomb self and exchange- 
correlation interactions of the localized d-orbitals of Sc have been corrected by the potential parameter of U. 
The EV cx -PW C o-GGA calculations have also provided good results for the structural and electronic features of 
ScN in ZB, WZ, and CsCI phases comparable with the theoretical data available in the literature. EVox-PWco- 
GGA and EVcx-PWco-GGA+U SIC schemes are considered to be the best ones among the others when the 
structural and electronic features of ScN are aimed to be calculated by the same exchange-correlation energy 
approximations. 
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1. Introduction 



ScN (Scandium nitride) has many potential applications due to its high mechanical strength [lfl, 
good thermal stability with a melting temperature above 2000° C [2j and high hardness of 21 GPa 
with respect to load deformation [l|. After the epitaxial growth of smooth and singly oriented ScN 
iy~j . films jlHg, ScN was also considered to be a potential semiconductor in electronic device applica- 

tions. X-Ray diffraction analysis indicated that, unlike the other nitride semiconductors, all ScN 
films crystallize into a single rock-salt (RS) phase (Bl) with a lattice constant of ~4.5 A [l|, [3h8|. 
The RS structure was also determined to be the most stable structure for ScN by the comparative 
theoretical works |9b11||. The lattice constant of ScN in RS phase has been calculated to be in 
the range of 4.42-4.651 A |9|. Illl4l5| by the first principles calculations within local density ap- 



proximation (LDA) [lg, ll7| and generaliz ed g radient approximation (GGA) [la . [l9j of exchange 
and correlation energies. In these works |9h14]- the bulk modulus of ScN has been reported in the 
range of 196-235 GPa with respect to the experimental value of 182=f40 [1|. In the literature the 
most stable RS phase with the cohesive energy of -13.69 [9j and -13.428 eV [10( was found to be 
followed by the wurtzite (WZ or B4), zinc-blende (ZB or B3), and CsCI (B2) type meta-stable 
structures of ScN with the ordered cohesive energies of -13.35 [9j, -13.03 [9fl, and -11.58 eV [9j 
(-11.34 eV [Uj), respectively. These cohesive energies correspond to the lattice constants of 3.49 
(c/a = 1.6, u = 0.38) @, 4.88 Hand 2.81 A |9] (2.79 A 0) for WZ, ZB, and CsCI phases of ScN, 
respectively. In a recent work [llj, the lattice constants have been calculated to be 3.45, 4.939, and 
2.926 A for ScN in WZ, ZB and CsClphases. The bulk modulus of WZ phase of ScN has been 
reported to be 156 [9j and 132.91 GPa [llj by GGA and LDA of exchange and correlation energies, 
respectively. In the same works, the bulk modulus of CsCI phase was calculated to be 170 [9j and 
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159.759 GPa [ll|. In another work [10| . the CsCl meta-stable phase of ScN has the bulk modulus 
of 178.60 GPa by GGA calculations. 

The first principles electronic band structure calculations within LDA [9j, [llj, [l2|, LLJ] and 



GGA [9j, [l0|) ll3| schemes indicated that RS-ScN was a semi-metal with an almost zero indirect band 
gap. However, ScN in RS phase was a semiconductor with the indirect gap of 2.4 [8j, 1.3=f0.3 [2GJ, 
and 0.9=f0.1 eV [7| measured by optical transmission, reflection and absorption experiments, re- 
spectively. The recent ab-initio calculations [12l [lj, [l5j, [2CJ, l21| using pseudopotential (PP) method 



within Hedin's [22j Green's functions quasiparticle corrections on LDA (LDA-G W ), exact ex- 
change of LDA [23| [OEPx(cLDA)], and Hedin's Green's functions quasiparticle corrections on 
exact exchange of LDA [0EPx(cLDA)-G o W o ], Linear Muffin Tin Orbital (LMTO) method within 
Bechstedt's [24] Green's functions quasiparticle corrections on LDA (LDA-GW) and full-potential- 
linearized augmented plane waves (FP-LAPW) method within GGA [19j and screened exchange 
of LDA (sx-LDA) [25| have all supported the semiconductor nature of the stable ScN b y g iving 
the indirect gap in the range of 0.54-1.70 eV at X symmetry point (E^~ x ). In [H, lTJ, SlHl , the 
optical direct gap of RS structure at X point (E X-X ) was calculated in the range of 1.98-2.90 eV 
with respect to the experimental direct absorptions in the range of 1.8-2.4 eV [3l.l4l.l7ll8ll20|. 

The recent DFT calculations within GGA scheme [9] have indicated that the WZ and ZB meta- 
stable structures of ScN were non-metallic with large indirect gaps of ~3 eV along M-r symmetry 
line (E^~ s ) and 2.3 eV at W symmetry point (E x_w ), respectively. The nonmetallic nature of 
ScN in ZB phase was also obtained by an indirect gap of 2.36 eV at W p oint (E x ~ w ) by LDA 
scheme [llj. These GGA and LDA calculations of ScN in ZB phase [9j, [llj, have given the direct 
gap of 2.4 and 2.42 eV at X symmetry point, respectively. In the literature, ScN in B2 phase was 
reported to be metallic by DFT calculations within GGA |9j, [lOj and LDA llj schemes. 



In the literature, although ScN has not been worked so far, hybrid FP-LAPW calculations 
within the framework of DFT and different exchange-correlation functionals have given accurate 
electronic features for nitride compounds and alloys 26J, [27 1 comparable with the corresponding 



measured ones due to the possible strong hybridization between the 2p orbitals of N and the 
corresponding cationic states. In addition, the orthogonalized norm-conserving pseudopotential 
(NCPP) method in which the plane wave basis functions are orthogonalized to core-like orbitals 
has been reported to be a very promising method in electronic band structure calculations of nitride 
alloys [28] for describing the experimental optical data, together with the FP-LAPW method within 
virtual crystal approximation [28]. 

In the present work, we have examined the structural and electronic properties of ScN in stable 
(RS) and meta-stable phases (CsCl, ZB, WZ) by DFT calculations mainly within two GGA schemes 
which have not been used before for ScN. We have aimed to introduce a comparative study for 
the structural and electronic features of ScN such as the lattice constant, bulk modulus, cohesive 
energy, energy gaps and the effective masses of electrons and holes. The present work has also 
comprised the electronic band structure of stable ScN corrected by an on-site Coulomb self- and 
exchange-correlation potential approximation (U SIC ) [29j. 

2. Method of calculations 

The present DFT calculations on the structural and electronic properties of ScN compound have 
been performed using FP-LAPW method implemented in WIEN2k code [30]. In the literature, the 
exchange-correlation energy of DFT has been defined by local density approximation (LDA) [16l | 
for the systems having uniform electron charge density. But for the systems of non-uniform charge 
density, the exchange-correlation energy of LDA has been corrected by gradient of the charge 
density within different generalized gradient approximations (GGA). In the present total energy 
and electronic band structure calculations of ScN in B1-B4 phases, four different approximations 
of exchange-correlation energies have been considered. In one of the approximations, exchange and 
correlation energies have been defined simply by LDA [16j , without regarding the homogeneity of 
the real charge density. In the second approximation, exchange and correlation energies of LDA have 
been corrected by GGA of Perdew-Burke-Ernzerhof (PBE) [19j. The generalized gradient functional 
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of Perdew-Burke-Ernzerhof [19] has retained correct features of LDA [16fl and satisfied only those 
which are energetically significant. In the third approximation, GGA of Engel-Vosko (EV) [3l| 
and GGA of Perdew and Wang (PW) [32j have been used to correct the exchan ge a nd correlation 
energies, respectively. Since the generalized gradient functional of Engel-Vosko |31| was designed 
to give a better exchange potential (V x ) only, the standard correlation potential of LDA [16j 
in the thir d ap proximation has been corrected by another functional, namely, GGA of Perdew 
and Wang [32|. The functional of Perdew and Wang [32] has incorporated some inhomogeneity 
effects while retaining many of the best features of the local density approximation. In the last 
approximation, the exchange energy of LDA [16j was corrected by GGA of Engel-Vosko [3l| , but the 
correlation energy was defined directly by LDA |l6| . The exchange-correlation energy approaches 
considered in this work have been labeled as LDA, PBE-GGA, EV cx -PW co -GGA and EV ex -GGA- 
LDA co . The acronyms have been produced either based on the key word of the approach (LDA) 
or on the name of the authors (PBE-GGA, EV ox -PW co -GGA, EV ex -GGA-LDA co ) who developed 
the corresponding exchange and correlation functionals. Here, the subscripts of exchange (ex) and 
correlation (co) functionals are exclusively used for the cases where the exchange and correlation 
functionals are different. It has been considered that EV ex -PW co -GGA and EV ex -GGA-LDA co 
schemes can provide significant improvement for the structural and electronic properties of ScN, 
respectively. In the literature, it was reported that LDA+U and GGA+U schemes can also improve 
the band gap energies of transition metal compounds and alloys by reproducing quite well the 
localized nature of the d-electrons (or /-electrons) [331435 1. Since the ab-initio calculations are 



difficult to perform, the strong correlations like in transition metal compounds and alloys are 
often based on a model Hamiltonian approach in which the important parameter of U improves 
the effective Coulomb interactions between the localized d-electrons. In the present work, the 
electronic band structure of RS-ScN calculated by EV cx -PW co -GGA has been improved by U 
method \M,\3$ introduced in WIEN2k code [30]. The present EV cx -PW co -GGA+U SIC scheme has 
rectified on-site Coulomb self- and exch ang e-correlation interactions of the localized d-orbitals of 
Sc by the potential parameter of U [29LI36J. In Anisimov et al's paper [29j, the meaning of the U 
parameter were defined as the cost Coulomb energy for the placement of two d-electrons on the same 
site. The Coulomb interaction was defined to be (1/2)U^ , nmj [37J for d-orbitals. Here, &j are 
d-orbital occupancies. The Coulomb interaction term included into the total energy functional of 
EV 0X -PW C o-GGA has given the orbital energies of ScN as ^= : C(EV el -PW 00 -GGA)+U(l/2 — m) [371 ]. 
The shifting of the corresponding orbital energies in EV cx -PW co -GGA+U SIC calculations gives a 
qualitative improvement for the energy gap of RS-ScN. The present Coulomb interaction parameter 
of U is calculated to be 4.08 eV in WIEN2k [30] to have the maximum approach to the measured 
indirect band gap of RS-ScN. 

In the present work, ScN has been studied in RS, CsCl, ZB, and WZ structures. The unit cells 
of RS, CsCl, and ZB consist of two basis atoms; Sc at (0, 0, 0) and N at (0.5a, 0.5a, 0.5a) in fee 
structure, Sc at (0, 0, 0) and N at (0.5a, 0.5a, 0.5a) in bec structure, and Sc at (0, 0, 0) and N at 
(0.25a, 0.25a, 0.25a) in fee structure, respectively, where a is the lattice constant parameter. The 
WZ structure with space group of F6^mc has four atoms in the unit cell; Sc atoms at (a/3, 2a/3, 0) 
and (2a/3, a/3, 0.5c), N atoms at (a/3, 2a/3, u*c) and (2a/3, a/3, (0.5 + w)*c), where a and c are 
the periods in x-y plane and along z direction, respectively. The z directional distance, u, is defined 
between the layers of Sc and N atoms. In FP-LAPW calculations, each unit cell is partitioned 
into non-overlapping muffin-tin spheres around the atomic sites. Basis functions are expanded in 
combinations of spherical harmonic functions inside the non-overlapping spheres. In the interstitial 
region, a plane wave basis is used and expansion is limited with a cutoff parameter, RmtKmax=7- 
Here, Rmt is the smallest radius of the sphere in the unit cell, if max is the magnitude of the 
largest K vector used in the plane wave expansion. The muffin-tin radius is adopted to be 1.8 and 
1.67 a.u. for Sc and N atoms, respectively. In the calculations, the electrons of Sc and N atoms in 
3s 2 3p 6 4s 2 3d 1 and 2s 2 2p 5 shells respectively, are treated as valence electrons by choosing a cutoff 
energy of -6.0 Ry. The core states are treated within the spherical part of the potential only and 
are assumed to have a spherically symmetric charge density totally confined inside the muffin-tin 
spheres. The expansion of spherical harmonic functions inside the muffin-tin spheres is truncated 
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at 1=10. The cutoff for Fourier expansion of the charge density and potential in the interstitial 
region is fixed to be Gmax = 16\/Ry. The FP-LAPW parameters presented in this work have been 
obtained after a few trials around their fixed values. 

The ScN in RS, CsCl and ZB structures have been optimized with respect to the volume of 
the unit cells by minimizing the total energy. The equilibrium lattice constants of ScN in RS, 
CsCl, and ZB phases are determined by fitting the total energies to the Murnaghan's equation 
of state [38j. The equilibrium structure of ScN in WZ phase that corresponds to the minimum 
total energy has been obtained by the application of both volume and geometry optimizations. 
The volume optimization used for all structures is provided with the energy criterion of 0.01 mRy. 
The optimum volume for WZ phase corresponds to the optimum c/a ratio and a value has been 
found by fitting the total energies to a quadratic function in a least square fitting method. The z 
directional distance u between the Sc and N layers in WZ phase has been obtained by geometry 
optimization at the optimum volume of the unit cell. The geometry optimization forces the atoms 
in the unit cell to move towards their equilibrium positions. In the geometry optimization, all forces 
on the atoms are converged to less than 1 mRy/a.u. The variation of total energy with respect 
to the volume and c/a ratio of WZ structure is plotted in figure [T] The present structural and 
electronic band calculations have been performed using 21x21x21 grids and correspond to 1000 k 
points sufficiently defined in the irreducible wedge of the Brillouin zone for ScN in B1-B4 phases. 



n = -3277.557040 Ry/unit cell 




80 Volume (a.u. ) 



Figure 1. (Color on-line) The relative total energy (per unit cell) versus volume and c/a within 
EV cx -PWco-GGA for WZ-ScN. 



In the present work, the cohesive energies (energy/atom-pair) of ScN structures (B1-B4) have 
been calculated by 



E. 



ScN 



ME n - ME Sc 



(2.1) 



Where, En and Es c are the values of the self-atomic energies of N and Sc atoms. Es c n is the 
minimum of the total energy per unit cell that corresponds to the equilibriunr structures of the 
ScN phases. Here, M defines the number of N and Sc atoms in the unit cell of the phases. The self- 
atomic energies of either atoms are calculated accurately in a fee super cell with a lattice constant 
of 25 a.u. 
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3. Results and discussion 

3.1. Structural properties 

The total energies (per unit cell) of RS, CsCl, WZ and ZB structures of ScN calculated within 
EV cx -PW co -GGA scheme are plotted as a function of the volume of the structures in figure [5] The 
curvature of variations is made clear by plotting the total energies (per unit cell) of the structures 
relative to the minimum total energy of the RS-ScN. The minimum of the total energies of the 
phases with E(RS-ScN)<E(WZ-ScN)<E(ZB-ScN)<E(CsCl-ScN) rank indicates that the RS is the 
most stable structure of ScN as it was reported in the other works |lLl3Hll| . The cohesive energies of 
ScN are calculated to be -12.40, -13.87, -14.06 and -14.48 eV for B2, B3, B4 type metastable and 
Bl type most stable structures, respectively. In the cohesive energy calculations [equation (2.1)], the 
atomic self-energies of Sc and N atoms are calculated to be -1528.609 and -109.136 Ry, respectively. 
The present cohesive energies are found to be different (~1 eV) from the ones |9j, [lOj defined by 
only the minimum total energies (per unit cell) of the structures. 
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Figure 2. The relative total energy (per unit cell) versus volume within EV cx -PW C o-GGA for 
ScN in RS, WZ, ZB, and CsCl structures. The inset figure displays the cross point on the 
variation of the total energies (per unit cell) of the RS-ScN (dotted line) and CsCl-ScN (solid 
line) structures as a function of volume. 



The equilibrium lattice constants, bulk moduli, and first order pressure derivative of the 
bulk moduli of ScN in B1-B4 phases calculated within LDA, PBE-GGA, EV cx -PW co -GGA, and 
EV cx -GGA-LDA co schemes are tabulated in table [TJ together with the available measured and cal- 
culated values of the other groups, for comparison. The equilibrium lattice constants of 4.510 and 
4.513 A obtained for RS-ScN structure by EV OX -PW CO -GGA and PBE-GGA calculations, respec- 
tively, are found to be very close to the measured values of 4.501 [lj and 4.5 A [3|-[8j|. On the other 
hand, the equilibrium lattice constant of RS-ScN obtained by the present EV cx -PW co -GGA and 
PBE-GGA calculations are found to be very close to the values calculated by FP-LAPW and ab- 
initio PP methods within PBE-GGA and PCVJPSF (Perdew-Chevary-Vosko-Jackson-Pederson- 
Singh-Fiolhais)-GGA J2^ schemes (HI H, E! • The present bulk moduli of 202.35 and 195.95 GPa 
found for RS-ScN within EV cx -PW co -GGA and PBE-GGA's, respectively, are the best ones pro- 
viding the measured value of 182=f40 GPa [lj. The present bulk moduli of EV ox -PW co -GGA 
and PBE-GGA's are also found to be in the range of 196-202 GPa defined by the correspond- 
ing results of ab-initio PP and FP-LAPW calculations within PBE-GGA and PCVJPSF-GGA 
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Table 1. The optimized lattice constant [a(A)j, bulk modulus [B(GPa)], and first-order pressure 
derivative of bulk modulus [B'] calculated by LDA (I), PBE-GGA (II), EV cx -PW co -GGA (III), 
and EVex-GGA-LDA co (IV) approaches for ScN in RS, CsCl, ZB, and WZ phases. 
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4.520i°, 4.651H 
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B 


159.12 
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142.35 
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a 
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3.494 
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3.499,3.45n 






c/a 


1.578 


1.582 


1.581 


1.581 


1.6 9 , 1.55211 




WZ-ScN 
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0.389 


0.389 


0.387 


0.38 9 ,0.377n 






B 
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131 


131 


86 


156 9 , 132.91H 






B 


4.18 


3.93 


3.43 


3.88 


2.16 9 ,3.29n 





schemes [9j, [lTJ, ll2Hl4 1 . The equilibrium lattice constants and bulk moduli of RS-ScN calculated 
by LDA and EV ex -GGA-LDA co 's are either underestimated or overestimated with respect to the 
measured ones given in table [T] The values of the first order pressure derivative of the bulk moduli 
(table [1} presented for RS-ScN are all close to each other and in agreement with the results of 
ab-initio PP calculations of LDA and PBE-GGA's [l3 |. As it is determined for ScN in RS phase, 
the structural features of ScN in B2-B4 phases calculated by EV ex -PW co -GGA and PBE-GGA's 
are very close to each other (table Q}. The present equilibrium lattice constants, bulk moduli, and 
first order pressure derivatives of the bulk moduli of ScN in B2-B4 phases could not be compared, 
because of the lack of the measured ones in the literature. But, they are (especially the results of 
LDA, PBE-GGA and EV BX -PW rn -GGA) comparable with the results of the other groups calculated 
within LDA [1, [Il| and PBE-GGA % [n| schemes. 

The total energies (per unit cell) of RS-ScN and CsCl-ScN relative to a reference total energy 
are plotted in the small frame of figure [2] for low unit cell volumes of the structures. The crossing 
of the total energy curves observed in the inset figure indicates the phase transition from RS to 
CsCl structure at high pressure. The enthalpy of the RS phase equals to that of the CsCl phase 



23701-6 



A comparative study for structural and electronic properties 




300 205 200 285 280 27-5 270 265 260 255 

Pressure (GPa) 

Figure 3. The difference of enthalpy between the CsCl and RS phases (AHcsCi-rs) versus 
pressure. 



at the cross point. Therefore, the necessary pressure providing the transition from RS to CsCl 
corresponds to the zero difference between the enthalpy of the structures. The enthalpy of the RS 
and CsCl phases is evaluated by Gibb's free energy 



H = E u 



PV-TS 



(3.1) 



at T=0 K. In figure[3j the difference between the enthalpy of the CsCl and RS structures (AiJcsCi-R.s) 
decreases to zero at the transition pressure of 285.2 GPa. The present transition pressure is found 
to be smaller than the ones (341 GPa [9|, 332.75 GPa [11|) calculated directly from the slope of 
the common tangent. 

3.2. Electronic properties 

The electronic band structures of ScN in RS, CsCl, ZB, and WZ phases have been calculated 
along the various symmetry lines within LDA, PBE-GGA, EV cx -PW co -GGA, and EV CX -GGA- 
LDA co schemes. The energy gaps of the RS-ScN correspond to different symmetry directions, and 
the points are tabulated in table[2j together with the available measured and calculated energy gaps 
of the other groups, for comparison. The electronic band structure calculations of RS-ScN phase 
have given the top of the valence band at T point. The bonding analysis given in [4CJ has shown 
that RS-A1N has three p-like bonding, three d-like antibonding t% s , and two d-like nonbonding 
e g bands formed by the hybridization of three- valence p states of N with the five d states of Sc. 
According to the present partial density of states (DOS) calculations plotted in figure [H the upper 
valence bands are formed mainly by the N p-orbitals with some mixture of Sc d-orbitals, while the 
conduction bands are predominantly originated from Sc d-t2 S states with some admixture of the Sc 
d-e g and N p-states. The minimum of the conduction band is mainly of Sc 3d character. The present 
identification for the conduction and valance bands of ScN is consistent with the partial density of 
states analysis given in [14J. The present LDA and PBE-GGA calculations have given the negative 
indirect gap at X point for RS-ScN. The negative or approximately zero indirect gap at X point 
was also reported in [9l4l4| for RS-ScN by the same approximations. The indirect energy gap of RS- 
ScN is 0.46 and 0.62 eV by the present calculations of EV cx -PW co -GGA and EV 0X -GGA-LDA co 's, 
respectively. Therefore, EV ex -PW co -GGA and EV ox -GGA-LDA co 's are found to be more accurate 
than PBE-GGA to produce a positive indirect gap for ScN in RS phase. These indirect gap values 
obtained by the present non-corrected EV ex -PW co -GGA and EV cx -GGA-LDA co calculations are 
in agreement with the value of 0.54 eV reported by non-corrected PP-GGA calculations [15]. In 
the literature, different correction approximations used in LDA and GGA schemes were found 
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Table 2. The calculated and experimental energy band gaps [E g (eV)j of RS-ScN. 
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Figure 4. (Color on-line) The total and partial DOS within EV cx -PW co -GGA for RS-ScN. 
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to be effective to produce a larger indirect gap close to the measured value of ~1 eV for RS- 



ScN [12|, UJ, [20l |2l|. The present EV ox -PW co -GGA+U SIC calculations have also improved the 



highest valence and lowest conduction band states of RS-ScN by giving an indirect gap of 0.9 eV. 
It is found that the corrected Coulomb interaction term by the U parameter of 4.08 eV pushes the 
N 2p bands down and Sc 3d-t2 g bands up in the present EV OX -PW CO -GGA+U SIC calculations to 
enlarge the indirect energy gap by an amount of 0.44 eV. Since the valence and conduction bands in 
RS-ScN structure are formed by the hybridization of d-orbitals of Sc and p-orbitals of N atoms [40|, 
the correction on the Coulomb self- and exchange-correlation interactions of the localized d-orbitals 
of Sc in EV cx -PW co -GGA+U SIC calculations has indirectly improved N 2p bands. In the same 
correction method, the direct band gap of RS phase at X point is calculated to be 1.82 eV with 
respect to the measured and calculated values in the ranges of 1.8-2.4 eV [3|, U, |7|, [20( and 1.98- 
2.90 eV [12l,llJ,l20j,l21|, respectively. The present T point direct energy gap of 3.28 eV for Bl phase is 
found to be comparable with the results of optical transmission [8] and reflection [20] measurements 
and 0EPx(cLDA)-G o W o [3] and LDA-G W Q calculations (table [S}. The present electronic 
band structure of RS-ScN within EV OX -PW CO -GGA scheme is plotted in figure[5l The improvement 
of the band structure due to EV OX -PW CO -GGA+U SIC calculations is also shown in the same plot. 



> 

•- 

o> 

S 




W K 



Figure 5. Electronic band structure of RS-ScN within EV C x-PW co -GGA. The dotted line shows 
the correction of the conduction band state within EV cx -PW C o-GGA+U calculations. 



In this work, the effective electron and hole masses (heavy and light) are calculated using 
the electronic band structure of RS-ScN obtained by EV ex -PW co -GGA+U SIC scheme (figure [SJ. 
The conduction band and valance band effective masses at the X and V points, respectively, are 
calculated by fitting a quadratic function to the corresponding band structure energies. The closely 
spaced k points in a very small range around the X point have yielded the electron effective mass 
of 0.22 m e . The electron effective mass in RS-ScN was determined between 0.1 and 0.2 m e by 
extrapolating the measured values of infrared reflectivity for high carrier concentrations to low 



4l|. Since the edge of the conduction band along A(r-X) and Z (W-X) directions approaches 



ones 

the X point with quiet different slopes (inset plot in figure [3]), the electron effective mass is also 
calculated along A and Z directions, separately. The electron effective masses are calculated to 
be 1.621 m e and 0.223 m e along A and Z directions around X point with respect to the ones 
reported in the range of 1.441-1.625 m e and 0.124-0.253 m e , respectively, by LDA 0,111, GGA 0, 
OEPx(cLDA) [14], and OEPx(cLDA)-G W [lj] schemes. In the present work, the effective mass 
for the heavy and light holes are calculated to be 0.93 and 0.19 m e at T point, respectively. 
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Figure 6. Electronic band structure of ZB-ScN 
within EV cx -PWco-GGA. 



Figure 7. Electronic band structure of WZ- 
ScN within EVox-PWco-GGA. 




r r x m r 

Figure 8. Electronic band structure of CsCl-ScN within EV ex -PW C o-GGA. 



The electronic band structures of ScN in B2-B4 metastable phases have been also calculated by 
LDA, PBE-GGA, EV OX -PW CO -GGA, and EV OX -GGA-LDA CO schemes. The band structures of the 
ZB, WZ, and CsCl phases calculated by EV ex -PW C o-GGA are plotted in figures [6HH The energy 
gaps of the ZB and WZ metastable structures correspond to different symmetry lines, and the 
points are tabulated in table [3] The calculations have yielded non-metallic band structures for 
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Table 3. The energy band gaps [E g (eV)j of ZB-ScN and WZ-ScN phases calculated by LDA (I), 
PBE-GGA (II), EV cx -PWco-GGA (III), and EV cx -GGA-LDA co (IV) approaches. 
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ZB-ScN and WZ-ScN with the top of the valence band at W and M symmetry points, respectively. 
The ZB-ScN structure is found to be a direct band gap material even if the direct (E^~ w ) and 
indirect (E^~ x ) band gaps are very close to each other. However, ZB-ScN was an indirect band 
gap (E x_w ) material in [9j, [llj. In reference [9j, the ZB-ScN structure was a direct band gap 
( E w-w-j ma t er j a i w hen the electronic band structures were calculated with the lattice constant of 
RS-ScN. The present direct band gap energies of PBE-GGA (2.44 eV) and LDA (2.30 eV) schemes 
for ZB structure are found to be comparable with the corresponding values of 2.4 and 2.36 eV 
given by similar approximations [9j, [llj. The present direct band gap of ZB phase at W point 
is enlarged to -3 eV by EV ox -PW co -GGA (3 eV) and EV ox -GGA-LDA co (2.82 eV) calculations. 
The indirect band gap of WZ phase is calculated to be ~3 eV along M-E direction by LDA and 
PBE-GGA schemes. The large indirect band gap of ScN in WZ structure was also reported in [9j. 
As it is found for ZB structure, the indirect band gap of WZ-ScN is enlarged by EV cx -PW co -GGA 
(3.29 eV) and EV ex -GGA-LDA co (3.08 eV)'s. According to all approximations considered in this 
work, ScN in CsCl phase has a metallic band structure (figure [5]) as it was reported earlier by 
LDA and GGA calculations [9j, [llj; the conduction and valance bands are observed to be mixed 
completely. Because of the lack of the measurements on metastable structures of ScN, the present 
results are compared only by the calculated ones. 

4. Summary and conclusion 

A comparative study by FP-LAPW calculations based on DFT within LDA, PBE-GGA, EVcx- 
PWco-GGA, and EV cx -GGA-LDA co schemes is introduced for the structural and electronic prop- 
erties of ScN in RS, CsCl, ZB, and WZ phases. According to all approximations used in this work, 
the RS phase is a stable ground state structure and makes a transition to CsCl phase at high 
transition pressure. It can be concluded that EV ex -PW C o-GGA is the best one among the others 
to provide accurate structural features and non-zero, positive indirect band gap of RS-ScN com- 
parable with the experimental results when the structural and electronic calculations are aimed to 
be calculated by the same exchange-correlation energy approximation. Although LDA and PBE- 
GGA's have calculated good structural features, they have yielded a small negative indirect band 
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gap for RS-ScN. On the other hand, EV cx -GGA-LDA co has roughly supplied the lattice constant 
and bulk modulus but it is found to be very accurate for the electronic features of RS-ScN. The 
indirect band gap of ScN in RS phase is enlarged to the corresponding measured value by EV ex - 
PW co -GGA+U calculations in which the Coulomb self- and exchange-correlation interactions 
of the localized d-orbitals of Sc have been corrected by the potential parameter of U. The present 
EV cx -PW co -GGA calculations have also provided good results for the structural and electronic 
features of ZB, WZ and CsCl phases comparable with the theoretical data reported in the lit- 
erature. Therefore, EV cx -PW co -GGA and EV cx -PW co -GGA+U SIC can be considered to be good 
exchange-correlation energy approximations for further works of ScN. 
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riopiBHsi/ibHe flocniflweiHhm CTpyKTypHux i e/ieKTpoHHux 

B/iaCTMBOCTeM MOHOKpHCTa/liHHOrO ScN 

P. MoraMMajiEl, C. KaTipHiorTiy^ 

riaiiecTWHCbKnii TexHiHHWw yHiBepcuTeT, KO/iefl>K npuK/iaflHMX HayK, BecTBeHK, riajiecTHHa 
BjiM3bKocxiflHw£i TexHiHHwii ymBepci/iTeT, cp i 3 i/i h h i/i w cpaKyjibTeT, 06530 AHKapa, TypeHHMHa 

ripeflCTaB/ieHO nopiBHajibHe flocniflweHHH 3a flonoMoroio FP-LAPW po3paxyHKiB, mo 6a3ytoTbOR Ha Teopii' 
cpyHKL4ioHa/iy rycTWHw (DFT) b paMKax cxeM LDA, PBE-GGA, EV cx -PW co -GGA, i EV cx -GGA-LDA co p/\n 
CTpyKTypHnx i ejieKTpoHHUx B/iacTi/iBOCTew ScN b cpa3ax RS, ZB, WZ i CsCI. BiflnoBiflHO flo Bcix Ha6/iw>KeHb, 
BMKOHaHnx b L4iti po6oTi, cpa3a RS e CTifiKoio CTpyKTypoio b ocHOBHOMy CTaHi i 3fli£icHK)e nepexifl y cpa3y 
CsCI npi/i BwcoKOMy TWCKy. B tow nac jrk cxeMW PBE-GGA i EV cx -PW C o-GGA 3a6e3nenyioTb Kpainj CTpy- 
KTypHi B/iacTWBOCTi TaKi ak piBHOBa>KHa nocTiKHa rpaTKH i 06'eMHi Mopym, nviine cxeMU EV cx -PW co -GGA Ta 
EV cx -GGA-LDA co 's flaioTb HeHy/ibOBy, no3WTHBHy HenpnMy eHepreTi/iHHy mi.ni/iHy p,nn RS-ScN, nopiBHajibHy 
3 eKcnepwMeHTajibHMMH. HenpaMa 30HHa mi^i/iHa ScN b cpa3i RS e 36mbiijeHa flo BiflnoBiflHoro BUMipnHoro 
3HaHeHHa 3a flonoMoroio EV cx -PW C o-GGA+U o6HwcneHb, b hki/ix KynomBCbKi B/iacHi i o6MiHHi Kope- 

jnsi L4iCi hi B3aeMOfliT jioKajii30BaHi/ix d-op6iTajie£n Sc 6y.ni/i nonpaB/iem 3a flonoMoroio napawieTpa noTeHL4ia.ny 
U. Po3paxyHKH EV cx -PW co -GGA TaKO>K npuBOflflTb flo flo6pwx pe3yjibTaTiB /\nn CTpyKTypHi/ix i ejieKTpoH- 
hi/ix xapaKTepucTWK ScN y cpa3ax ZB, WZ i CsCI, hklho nopiBHioBaTi/i 3 TeopeTi/iHHUMW flaHMMW, HaflBHUMi/i 
b jiiTepaTypi. BBaxoeTbcsi, mo cxeMi/i EV cx -PW co -GGA Ta EV cx -PW co -GGA-(-U e HawKpamnMH cepefl 
iHwwx y BunaflKy, ko/ih o6HncntotoTbc:R CTpyKTypHi Ta ejieKTpoHHi xapaKTepwcTHKH ScN b pawiKax tux we 
Ha6.ni/i>KeHb fl/ifl eHepriT o6MiHHoT Kopejiai^iT. 

KyiKJHOBi c/iOBa: ScN, FP-LAPW, DFT, CTpyKTypHi B/iacTMBOcri, e/ieKTpoHHi BnacTWBOCTi 
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